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Abstract 

The probability density function (PDF) of accelerations in turbulence is derived 
analytically with the help of the multifractal analysis based on generalized entropy, 
i.e., the Tsallis or the Renyi entropy. It is shown that the derived PDF explains quite 
well the one obtained by Bodenschatz et al. in the measurement of fluid particle 
accelerations in the Lagrangian frame at Rx = 690, and the one by Gotoh et al. in 
the DNS with the mesh size 1024^ at Rx = 380. 
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1 Introduction 



The multifractal analysis of turbulence [1,2,3,4,5,6,7,8,9,10] is a unified self- 
consistent approach for the systems with large deviations, which has been con- 
structed based on the Tsallis-type distribution function [11,12] that provides 
an extremum of the extensive Renyi [13] or the non-extensive Tsallis entropy 
[11,12,14] under appropriate constraints. The analysis rests on the scale in- 
variance of the Navier-Stokes equation for high Reynolds number, and on the 
assumptions that the singularities due to the invariance distribute themselves 
multifractally in physical space. The multifractal analysis belongs to the line 
of study based on a kind of ensemble theoretical approaches that, starting 
from K41 [15], continues with the log-normal model [16,17,18], the /3-model 
[19], the p-model [20,21], the 3D binomial Cantor set model [22] and so on. 
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After a rather preliminary investigation of the p-model [1], we developed fur- 
ther to derive the analytical expression for the scaling exponents of velocity 
structure function [2,3,4,5], and to determine the probability density function 
(PDF) of velocity fluctuations [5,6,7,8], of velocity derivative [9] and of fluid 
particle accelerations [10] by a self-consistent statistical mechanical approach. 
It has been shown [5] that the multifractal analysis derives the log-normal 
model [16,17,18] when one starts with the Boltzmann-Gibbs entropy. 

In this paper, we derive the formula for the PDF of the accelerations of a fluid 
particle in fully developed turbulence by means of the multifractal analysis, 
and will analyze two experiments. One is the PDF of accelerations at Rx = 
690 of the Taylor microscale Reynolds number obtained in the Lagrangian 
measurement of particle accelerations that was realized by Bodenschatz et 
al. [23,24,25] by raising dramatically the spatial and temporal measurement 
resolutions with the help of the sihcon strip detectors. The other is the PDF 
of accelerations at Rx = 380 [26] that was extracted by Gotoh from the DNS 
of the size 1024^ [27] which may be the largest mesh size available at present. 

For high Reynolds number Re » 1, or for the situation where effects of the 
kinematic viscosity u can be neglected compared with those of the turbulent 
viscosity, the Navier-Stokes equation, 

du/dt + {u ■ V)u = -V (p/p) + uV^u, (1) 

of an incompressible fluid is invariant under the scale transformation [28,29,21] 
r^Af, u^X''/^ t^X^-'^/H, (p/p) ^ A'"/3 (p/p) (2) 

where the exponent a is an arbitrary real quantity. The quantities p and p 
represent, respectively, mass density and pressure. The acceleration a of a fluid 
particle is given by the substantive time derivative of the velocity: 

a = du/dt + {u ■ V)u. (3) 

The Reynolds number Re of the system is given by 

Re = 6uaiji^ = {i-JviY'^ (4) 

with the Kolmogorov scale 77 = [u'^ /eY^^ [15] where e is the energy input rate 
at the energy-input scale £in. The velocity fluctuation 5u\n is defined by putting 
= (-in in 5un = \u{» + £n) — «(•) | where w is a component of velocity field u, 
and in is a distance between two observing points. We are measuring distance 
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by the discrete units 



(5) 



with 5n = 2~" [n = 0, 1, 2, • • •) and ^o being some reference length. The non- 
negative integer n represents the multifractal depth. However, we will treat it 
as positive real number in the analysis of experiments. 

It may be worthwhile to note here that, within the energy cascade model, the 
diameter in of eddies should be defined by 

4 = SJin. (6) 



Therefore, if one identifies the distance with the diameter 4 of the nth 
eddies within the energy cascade model, one has the relation 

ri = n-log2(4Ai„) (7) 



between the number n of the multifractal steps and the number n of the steps 
within the energy cascade. We see that £in = io ^ it should be. 

Introducing the pressure (divided by the mass density) difference 

Spn^\p/P{*+Q-P/Pi*)\ (8) 

between two points separated by the distance in, the accelerations may be 
estimated by 

lal = lim a„ (9) 

n— >cxD 

where we introduced the acceleration 

an = Spn/in (10) 

belonging to the multifractal depth n. The acceleration becomes singular for 
a < 1.5, i.e., 

lim a, = )im Spn/in ~ lim/(f"/=^)-^ ^ oo (11) 

which can be seen with the relation 

Spn/Spo = iL/iof"^'. (12) 
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The values of exponent a specify the degree of singularity. 



2 Multifractal spectrum 



The multifractal analysis rests on the multifractal distribution of a. The prob- 
ability P^"\a)da to find, at a point in physical space, a singularity labeled by 
an exponent in the range a ~ a + da is given by [2,3,4,5] 

p(-\a) = \l-{a- aor/{Aa)T^^'~'^ (13) 



with an appropriate partition function and 

{Aaf^2X/[{l-q)]n2]. (14) 

The range of a is ccmin < cc < a^ax with 

ttmin = Q!0 - Aa, Ctmax = "O + Aa. (15) 

It may be worthwhile to put here its brief derivation in order to make the 
paper self-contained. The distribution function (13) is derived by taking an 
extremum of the generalized entropy, the Renyi entropy [13] 

Sf[P^^\a)] = (1 - q)-^ In J daP^^\ay (16) 

or the Tsallis entropy [11,12,14] 

Sj[P^'\a)] = (1 - q)-' (I da P^'\aY - l) , (17) 

under the two constraints, i.e., the normalization of distribution function: 
J daP^^\a) ^ const., (18) 

and the g-variance being kept constant as a known quantity: 

{ay {a — ao f] daP^^\ay. (19) 
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Here, we assTime that the distribution function at the nth multifractal depth 
has the structure 



pin) 

(a) oc [P^^\a)\^. (20) 

This is consistent with the relation [21,5] 

P(")(«) oc (5^^^°) (21) 



that is a manifestation of scale invariance and reveals how densely each singu- 
larity, labeled by a, fills physical space. In the present model, the multifractal 
spectrum f[a) is given by [2,3,4,5] 

/(a) = 1 + (1 - q)-' log2[l - (a - ao)V(A«)']- (22) 



In spite of the different characteristics of these entropies, i.e., extensive and 
non-extensive, the distribution functions giving their extremum have the com- 
mon structure (13). ^ 

The dependence of the parameters ao, X and q on the intermittency expo- 
nent // is determined, self-consistently, with the help of the three independent 
equations, i.e., the energy conservation: 

(en) = e, (23) 



the definition of the intermittency exponent ji: 

(el) = e%^, (24) 

and the scaling relation ^ : 

l/{l-q) ^l/a--l/a+ (25) 

with a± satisfying f{a±) — 0. The average (• • •) is taken with P^^^ia). The 
energy-transfer rate e„ represents the rate of transfer of energy per unit mass 
from eddies of size In = in (the nth step within the energy cascade model) to 
those of size = ^7^^ (the n + 1th step within the energy cascade model) . 



^ Within the present formulation, the decision cannot be pronounced which of the 
entropies is underlying the system of turbulence. 

^ The scahng relation is a generalization of the one derived first in [30,31] to the 
case where the multifractal spectrum has negative values. 
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It should be noted that the average (e„) is taken just for the eddies having 
the size £„. It is not the average within the spatial region of the diameter 

For the region 0.13 < // < 0.40 where the value of pL is usually observed, the 
three self-consistent equations are solved to give the approximate equations [8] : 

ao = 0.9989 + 0.5814/^, X = -2.848 x 10"^ + 1.198/x (26) 

q = -1.507 + 20.58/x - 97.11f/ + 260.4//^ - 365.4/x^ + 208.3//^ (27) 



3 Scaling exponent of velocity structure function 



Let us derive first the probability A*^")(y„)(iy„ to find the scaled pressure fluc- 
tuations 

\yn\= 5pn/5po (28) 



in the range i/n yn + dUn in the form 

A^"^ {yn)dyn = ^^t^ {yn)dyn + AA(") {yn)dyn (29) 

with the normalization 

%nA(")(l/n) = 1- (30) 



Here, we assumed that it has two contributions whose origins are independent 
with each other. The first term represents the contribution by the singular part 
of accelerations stemmed from the multifractal distribution of the singularities 
in physical space. This is given by 

A^s\\yn\)dy„<xP^^\a)da (31) 



with the transformation of the variables, \yn\ = Whereas the second 

term AA^")(y„)dy„ represents the contribution from the dissipative term in 
the Navier-Stokes equation, and/or the one from the errors in measurements. 
The dissipative term has been discarded in the above investigation since it 
violates the invariance under the scale transformation. The contribution of 
the second term provides a correction to the first one. Note that the propor- 
tionality coefficient in (31) determines the portion of the contribution among 
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these two independent origins. ^ 

The mth moments of the pressure fluctuations, {{\yn\™')), are given by 

OD 

{{\yX')) ^ I dyn\ynr^^'\yn) = ^it^ + (1 - 2#) a^^ 5^- (32) 



where 



27H = I |y„rAAH(y„) (33) 
—00 

«3, = {2/[^(l + ^)]r/^ (34) 

with 

Q= 1 + 2^2(1 -g)X In 2. (35) 
The quantity 

= aom/3 - 2XmV[9(l + C%)\ - [1 - log2(l + Ci/;3)]/(l - q) (36) 

is the scahng exponent of the mth order velocity structure function which 
explains successfully the experimental results [2,3,4,5]. Note that the formula 
is independent of n. 

4 PDF of the fluid pcirticle accelerations 

We now derive the PDF, A('*)(a;„), defined by the relation 

A(")(a;„)du;„ = A(")(y„)dy„ (37) 

with the acceleration lo^ normalized by its deviation, i.e., 

W\ = a„/((a^))V2 = hn^mY'^ = Vnliy'S'' = u^Jl-l'-'^l^ (38) 



^ Needless to say that each term in (29) is a multiple of two PDFs, i.e., the PDF for 
one of the two independent origins to realize and the conditional PDF for a value 
y„ in the range ?/n ~ 2/n + <^2/n to come out. This is of course in a generalized sense 
in which the second correction term may weaken the first singular contribution. 
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with 



= [2lt^5-^* + (1 - 27rO«4, 



(39) 



It is reasonable to imagine that the origin of intermittent rare events is at- 
tributed to the singular term in (29). We then have, for the tail part uj^^ < 

\, , \ ^ , .max 
l^nl S 1 



(n) 



I- q (3 In \Un/^^nfi 



2-1 n/(l-g) 



n 



8X1 In (5. 



with 



^ ^ ^ ^2ao/3-C4/2 max ^ - ^2a^i„/3-C4/2 



(40) 



(41) 



and 



A(") = 3 (l - 27$-)) / (4u;„^27rX|ln5. 



(42) 



On the other hand, for smaller accelerations, the contribution to the PDF 
comes from both the singularity and thermal fluctuations or measurement er- 
ror. We assume that this part of the PDF is described by Tsallis-type function 
with a new parameter g', i.e., for the center part \u}n\ < u;^. 



k^^\uJn)dUn = \kt\yn) + AAW(y„)l dy, 



l-q' 



1 + \no.') 



^1 _ 1 



1/(1-9') 



dWn- 



(43) 



This specific form of the Tsallis function is determined by the condition that 
the two PDFs (40) and (43) should have the same value and the same slope 
at UJl^ which is defined by 



(44) 



with being the smaller solution of 
C4/2 - 2q;/3 + 1 - /(a) = 0. 



(45) 



It is the point at which A(")(c<;^) has the least n-dependence for large n. 
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With the help of the relation (43) , we obtain AA^") , and have the formula 
to evaluate in the form 



(46) 



where 



3 ^2(m+l)Qt/3_^4/2 1 



2^2nX\lndn\ { 



dz z"' 



X 



1 - 



l-q' 



1 + ^/'(«^) 



(.^ - 1) 



1/(1-5') 



2J27rX|ln5, 



^ _ l-q (3 ln|z/4|)' 
n 8X|ln(5„| 



n/{l-q) 



(47) 



(48) 



with 



^ _ , , /, ,t — x2{amax-Q;''')/3 _ , , /, ,t _ x2{Qo-at)/3 



(49) 



Now, the PDF of fluid particle accelerations, given by (40) and (43), is com- 
pletely determined by three parameters, i.e., the intermittency exponent /x, 
the multifractal depth n which gives a characteristic length and q' which 
appears in the Tsallis-type PDF at the center part. The intermittency expo- 
nent /i is determined by analyzing the measured scaling exponent Cm of the 
velocity structure function with the formula (36) as mentioned before. 

As for the determination of the value n, the flatness of the PDF of fluid particle 

accelerations can be a good candidate [25]. Actually, the present multifractal 
analysis provides us with the analytical formula for the flatness of the PDF in 
the form 



27^) 



as 



+ 



as 



^C8-2C4 

:i- 27^)^1 " ■ 



(50) 
(51) 



With this formula, the experimental value of the flatness F^"-'> gives the value 
of the multifractal step n. In deriving the approximate formula (51), we used 
the fact that the first terms both in the denominator and numerator of the last 
formula in (50) are two or three orders in magnitude smaller than the second 
terms. Note that the contribution to the flatness is mainly come from the tail 
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Fig. 1. Dependence of the flatness on the multifractal step n. The thin and thick 
hnes represent, respectively, the dependence for the parameters used in the study 
of the experiment by Bodenschatz et al. (Fig. 2 below), and of the DNS by Gotoh 
et al. (Fig. 3 below). 

part of the PDF, and that there is almost no contribution from its center 
part (see Fig. 4 below). Therefore, at this stage, we can determine almost 
the final shape and the magnitude of the tail part given by (40). It means 
that the shape of the center part, which is controlled by q', rarely affect the 
tail part, and that the contribution of the center part to the normalization of 
PDF is almost independent of q\ The n-dependence of the flatness is given 
in Fig. 1 both for the cases corresponding to the experiments conducted by 
Bodenschatz et al. with the approximate formula 

F^") = -6.872 + 5.216n - 9.969 x lO^^n^ 

+1.075 X 10"V - 5.296 x 10"V + 1.241 x 10"V, (52) 

and to the DNS performed by Gotoh et al. with 

F^H = -7.592 + 5.831n - l.lMn^ 

+1.203 x 10- V - 5.925 x lO^V + 1.389 x lO'V. (53) 

The slight difference between the two formulae comes from the difference of 
the values q' (see Fig. 2 and Fig. 3 below). 

Within the present stage of the multifractal analysis, the value q' is determined 

by adjusting the center part with the experimental data. However, since the 
center part of the PDF takes care only smaller accelerations compared with 
its deviation, we can expect that the dynamical analysis such as done by 
Beck [32,33] may be effective. The dynamical study for the center part will be 
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Fig. 2. PDF of accelerations plotted on (a) log and (b) linear scale. Comparison 
between the experimentally measured PDF of fluid particle accelerations by Boden- 
schatz et al. at Rx = 690 (Re = 31 400) and the present theoretical PDF A(")(a;„). 
Open squares axe the experimental data points on the left hand side of the PDF, 
whereas open circles are those on the right hand side. Solid lines represent the curves 
given by the present theory (40) and (43) with = 0.240 {q = 0.391), n = 17.1 and 
q' = 1.45. 

reported elsewhere in connection with the multifractal characteristics of the 
system which determines the tail part representing large deviations. 



5 Analysis of experiments 

We analyze, with the formula (40) and (43), the experimental PDF at Rx — 
690 measured by Bodenschatz et al. [25] in Fig. 2 on (a) log and (b) linear 
scale. We determined the value n — 17.1 for this experiment by substituting 
into its definition (5), i.e., n = log2(£o/0) the reported value of the integral 
length scale 0.071 m for £q, and of the spatial measurement resolution 0.5 fim 
for in [23,24]. The intermittency exponent ii = 0.240 is extracted by the 
method of least squares with respect to the logarithm of PDFs as the best fit 
of our theoretical formulae with n — 17.1 to the observed values of the PDF 
[23,24,25]. We discarded those points whose PDF values arc less than ~ 10~^ 
since they scatter largely on log scale. Substituting the extracted value of fi 
into the self-consistent equations, wc have the values of parameters: q = 0.391, 
ao = 1.138 and X = 0.285. With these values, other quantities are determined, 
e.g., Aa — 1.161, a+ — ao — ao — a- — 0.6815, a^ — 1.005, o;^ = 0.605 and 
^max ^ 2040. The best fit at the center part is given with = 1.45.'' The 

^ It is remarkable that this value of q' is close to q' = 1.5 proposed by Beck [32,33], 
although it is claimed that his formula does not provide a good explanation of 
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Fig. 3. PDF of accelerations plotted on (a) log and (b) linear scale. Comparison 
between the PDF of fluid particle accelerations measured in the DNS by Gotoh 
et al. at Rx = 380 and the present theoretical PDF A(")(a;„). Closed circles are 
the DNS data points both on the left and right hand sides of the PDF. Solid lines 
represent the curves given by the present theory (40) and (43) with /i = 0.240 
{q = 0.391), n = 17.5 and q' = 1.70. 

flatness of the PDF of accelerations has the value F^^^ = 56.9 which is within 
the reported value of the flatness 55 ± 4 [25]. This results tells us that the 
formula (50) of the flatness can be used to derive the value of /i, accurately, 
when one knows the value of n, and vice versa (see Fig. 1 and also Fig. 4 
below) . 

Analysis of the PDF of accelerations extracted out from the DNS data ob- 
tained by Gotoh et al. [27] at R\ = 380 is shown in Fig. 3 on (a) log and 
(b) hnear scale. The value fj, — 0.240 for this DNS has been determined by 
the analysis of the experimental scaling exponent Cm of longitudinal velocity 
structure function with our theoretical formula (36) [8,9]. Substituting this 
value of /i into the self-consistent equations, we have the values of parameters: 
q = 0.391, ao = 1.138, X = 0.285, Aa = 1.161 and a+-ao = ao-a. = 0.6815 
which are the same as those derived in the analysis of the experiment by Bo- 
denschatz, since both systems have a common value fi — 0.240 for the inter- 
mittency exponent. The value n = 17.5 is extracted by the method of least 
squares with respect to the logarithm of PDFs as the best fit of our theoret- 
ical formulae with the derived parameters given above to the observed data 
of the PDF [27]. Note that = 1.005, cul = 0.622 and = 2534. The 
best fit at the center part is given with q' — 1.7. The fiatness of the PDF 
is F^"-^ = 70.0. The characteristic distance r = in for n = 17.5 reduces to 
r/rj = 7.91 with 1] ^ 0.258 x 10^^ [27], which is about three times longer 
than the length Ar/r] = 27r/(1024r/) = 2.38 of the mesh of Gotoh's DNS. In 



experimental PDF for larger accelerations [25]. 
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deriving the characteristic distance, we used the second formula in ^ 

n = -1.050 log2(r/r7) + 16.74 (4 < r) (54) 
n = -2.540 log2(r/77) + 25.08 (r < 4). (55) 

Here, ic/v = 48.26 is the crossover length of two scaling regions, and is close 
to the Taylor microscale X/r] = 38.33. With the help of the formula (7), we 
have n — 11.8. 

In Fig. 4, we put the lines representing the integrand u^A^^^ (un) of the flatness 
and the corresponding experimental data (a) by Bodenschatz et al. [25], 
and (b) by Gotoh et al. [27]. The agreements are remarkable. 



6 Comment on the energy-input scale 

There is no room to incorporate, automatically, into the present multifractal 
analysis the energy input scale and the "system size" i^. The former is nec- 
essary to determine the number of steps n in the energy cascade model. Once 
£in is determined by investigating the structure of experimental apparatuses, 
the relation between n and the multifractal step n is given by (7). Since main 
part of the multifractal analysis rests on the scale invariance, the size of the 
system under consideration is assumed to be infinite, and therefore, the value 
of the reference length io introduced in (5) is determined only through the 
analysis of experimental data. 

Actually, for example, the empirical equation [6] 

n = -1.019 log2r/r/ + 0.901 log2 Re (56) 



^ The empirical formulae (54) and (55) were extracted through the multifractal 
analysis for the PDFs of velocity fluctuations measured at two points separated by 
r [8]. The corresponding numbers n and n are given by (r/jy, n, n) = (2.38, 21.5, 
14.5), (4.76, 20.0, 13.0), (9.52, 16.8, 9.81), (19.0, 14.0, 7.01), (38.1, 11.8, 4.81), (76.2, 
10.1, 3.11), (152, 9.30, 2.31), (305, 8.10, 1.11), (609, 7.00, 0.01), (1220, 6.00, -0.99). It 
seems that there arc two scaling regions constituting the inertial range. The formula 
(55) indicates that in the lower scaling region (r < 4) eddies break up, effectively, 
into 1.33 pieces in contrast to the upper scaling region (4 < r) where eddies break 
up into two pieces as seen in (54) [8]. The existence of the two scaling regions can 
be an artifact of DNS due to a shortage of calculation time, since energy cascading 
process is not so effective for eddies whose sizes are smaller than ic or the Taylor 
microscale. Note that Gotoh et al. claimed that the inertial range is restricted only 
to the region around n = 3.11, 2.31 and 1.11 [27]. 
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Fig. 4. Comparison of the theoretical curves (sohd lines) for the integrand 
ujf^M^\uJn) of the flatness F^"'^ with the corresponding (a) experimental quantity 
by Bodenschatz et al. (open squares for the data on the left hand side, and open 
circles on the right hand side), and (b) DNS quantity by Gotoh et al. (closed circles 
both for the data on the left and right hand) in loglog scale. The parameters in the 
theoretical PDFs for (a) and (b) are the same as those given in Fig. 2 and Fig. 3, 
respectively. 

extracted from the experimental PDFs for velocity fluctuations^ by Lewis 
and Swinney [34] gives £o ~ 877 cm. In their analysis, the Reynolds number 

Re = 540 000 is estimated with 4^ = 27r x 19.00 cm ^ 119.32 cm and the Kol- 
mogorov scale rj ^ 0.006 cm. Therefore, the reference length is large compared 
with the energy-input scale. 

For Gotoh's DNS, the empirical equation (54), extracted by the analysis of 
PDFs of velocity fluctuations, gives io/r] ^ 63 000 which is larger than the 
energy-input scale iin/i] = 7r/(kr]) ^ 497 with the wavenumber k = 6^^^ of 
forcing where we took 209 grid point spacings for the measure of the energy- 
input scale [35]. We see again that the reference length is longer than the 
energy-input scale. The Reynolds number of the DNS is now estimated by 
(4) as Re = 3 937. A unifled study of Gotoh's DNS [27] by means of the 
multifractal analysis with the Tsallis-type PDF for the center part will be given 
elsewhere, in which it is shown that the PDFs for the velocity fluctuations, for 



The velocity fluctuations are measured at two separated points whose distance is 
r. The corresponding numbers n and n are given by [10] {r/rj, n, n) = (11.6, 14, 
10.7), (23.1, 13, 9.7), (46.2, 11, 8.7), (92.5, 10, 7.7), (208, 9.0, 6.6), (399, 8.0, 5.6), 
(830, 7.5, 4.6), (1440, 7.0, 3.8). The relation (56) is extracted with these values [6]. 
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the velocity derivatives and for the fliiid particle accelerations in addition to 
the scale exponents of velocity structure function provide us with consistent 
results. 

In the analysis of Bodenschatz's experiment, the identification of io with the 
integral length scale 7.1 cm gives us reasonable value of n = 17.1 in the 
sense that with this value the formula (50) of the flatness provides us with 
Ij, = 0.240 for the intermittency exponent. This value turns out to be the same 
as the one observed by Gotoh et al. in their DNS. Since Bodenschatz assigned 
Re = 31 400 with the formula (4), rj — 30.3 //m gives £in = 0.071 m, as it 
should be. Therefore, in this case, we see that the energy-input scale is equal 
to the integral length scale, i.e., ii^ = £o, and that n = n = 17.1. 



7 Discussions 

A comparison of the PDFs of accelerations are put in Fig. 5 on (a) log and (b) 
linear scale. The thin and the thick lines are those theoretical PDFs in Fig. 2 
and in Fig. 3, respectively. The analysis of the data by Bodenschatz for the 
accelerations in terms of the log-normal model was performed just the same 
way as the multifractal analysis to get // = 0.260 and n — 15.5. The resulting 
PDF is given by dashed line which deviates slightly upward for larger values 
of ujn compared with the thin line. This deviation can be understood from 
the scaling exponents (m of the velocity structure function within the log- 
normal model. It becomes negative for larger values of m. The dotted hne is 
the empirical PDF given by Bodenschatz et al. [25] which deviate downward 
for larger values of cj„. The dotted-dashed line is the PDF given by Beck [36] 
for the studies both of the PDFs by Bodenschatz et al. and of Gotoh et al.. 
The lines other than the present theoretical PDFs in Fig. 5 give rather poor 
explanations at the center part. ^ 

From the above analyses of two experiments, we reveal that there are two 
mechanisms contributing to the PDF of the accelerations, i.e., one is for the 
tail part, and the other for the center part. The structure of the PDF A*^"^(a;„) 
for the tail part, a;^ < \cOn\, is given by (40) that represents the intermit- 
tent large deviations which is a manifestation of the multifractal distribution 
of singularities in physical space due to the scale invariancc of the Navier- 
Stokes equation for large Reynolds number. The experiment conducted by 
Bodenschatz et al. [23,24] visualized the singularities in physical space by 
tracing the "fluid particle" in turbulence. The speciflc form (40) comes from 
the distribution function for the singularity exponent a that is represented 

^ Beck introduced a somewhat artificial cutoff in his formula in order to fix its 
defect at the center part of the PDF [37]. 
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Fig. 5. Comparison of PDFs of accelerations on (a) log and (b) linear scale. The 
thin and thick lines are, respectively, the present theoretical PDFs in Fig. 2 and 
Fig. 3. The dashed line represents the PDF within the log-normal model, and the 
dotted line docs the empirical PDF given by Bodenschatz et al.. The dotted-dashed 
line is the PDF given by Beck [36]. 

by the Tsallis-type distribution function (13) with the parameter q which is 
determined by the observed value of the intermittency exponent. The flat- 
ness of the PDF mainly provides us with the information of this tail part. 
The structure of the PDF for the center part, |ci;n| < u]^, is given by (43) 
that represents small deviations violating the scale invariance due to thermal 
fluctuations and/or observation error. The PDF for this part is assumed to 
be given by the Tsallis-type distribution function for acceleration itself with 
the parameter q' . In this paper, the value of q' are determined with the help 
of the experimentally observed PDF at the center part, giving q' = 1.45 and 
q' = 1.7 which is close to 1.5 proposed by Beck [32,33]. The value of q' should 
be determined by investigating the dynamics of thermal fluctuations and/or 
observation error on the multifractal support, which may have non-additive 
character. This is one of the attractive future problems. Note that we already 
saw a dependence of ln(l/(g' — 1)) on \n{r/ri) through the study of the PDFs 
of velocity fluctuations, a detailed of which will be reported elsewhere. The 
tail part and the center part are separated at a;^ ~ 0.6. This gives ^ 1 
that satisfles the condition a < 1.5 in which the singularity appears in fluid 
particle accelerations. 

We saw that the PDFs derived within the multifractal analysis seem to be 
sensitive to the characteristic lengths such as the distance of two measuring 
points, the space resolution in measurement and the mesh size of DNS. We 
also knew that the multifractal distribution of singularities in physical space, 
on which the present analysis rests, is robust enough to allow us to apply the 
multifractal analysis to the ranges outside of the inertial range. How to put the 
information of energy-input scale into the multifractal analysis is one of the 



16 



important future problems [38]. It may be resolved when one succeeds to reveal 
the dynamical foundation underlying the basis of the multifractal analysis, 
starting an investigation by the stochastic Navier-Stokes equation with the 
energy input term. If one could put the stochastic Navier-Stokes equation 
under the influence of white Gaussian noise, describing thermal fluctuation 
related to kinematic viscosity, into a linealizcd stochastic equation with a 
renormalizcd turbulent viscosity [39], the relevant stochastic process should be 
the one related to the PDF of accelerations derived in this paper, which may be 
named Renyi or Tsallis process. The success of the multifractal analysis in the 
studies of experimental PDFs may indicate that the multifractal distribution 
of singularities in physical space are robust against the addition of the energy- 
input term. 

The application of the present multifractal analysis to the micrometeorological 
study of atmospheric turbulence is one of the attractive future problems. The 
tail part is specified by the intermittency exponent, and the center part may 
give us the information, for example, about the thermodynamical structure 
of the canopy in Amazon forest [40]. The application to the vortex tangle 
[41] in the superfluids ^He and ^He is another examples of the exciting future 
problems. In low temperature, the tail part of PDF comes from the singularity 
in the supcrfluid component within the two fluid model, and is determined by 
the intermittency exponent. Whereas, the value of q' for the center part of 
PDF may be determined by the dynamical structure of quantized vortices. 
The multifractal analysis may also open new aspect for the systems with large 
deviation found in a large variety of areas (see, for example, the web site in 
[12]), and its apphcations to these areas will be reported elsewhere. 
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